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Abstract 

We consider an initial-boundary value problem for a fully nonlinear coupled parabolic sys- 
tem with nonlinear boundary conditions modelling hygro-thermal behavior of concrete at high 
temperatures. We prove a global existence of a weak solution to this system on an arbitrary 
time interval. The main result is proved by an approximation procedure. This consists in 
proving the existence of solutions to mollified problems using the Leray-Schauder theorem, for 
which a priori estimates are obtained. The limit then provides a weak solution for the original 
problem. A practical example illustrates a performance of the model for a problem of a concrete 
segment exposed to transient heating according to three different fire scenarios. Here, the focus 
is on the short-term pore pressure build up, which can lead to explosive spalling of concrete 
and catastrophic failures of concrete structures. 

1 Introduction 

The hygro-thermal behavior of concrete plays a crucial role in the assessment of the rehability and 
lifetime of concrete structures. The impact of heat and mass transfer processes become particularly 
evident at high temperatures, where the increased pressure in pores, large temperature gradients 
and temperature induced-creep may lead to catastrophic service failures. Examples of such severe 
situations include hypothetical nuclear reactor accidents [30], fires in tunnels and tall buildings 
|31) or simulation of airfield concrete pavements |21) . Since high-temperature experiments are very 
expensive, predictive modeling of humidity migration, pore pressure development and damage dis- 
tribution can result in significant economic savings. The first mathematical models of concrete 
exposed to temperatures exceeding 100°C were formulated by Bazant and Thonguthai in [S]. Since 
then, a considerable effort has been invested into detailed numerical simulations of concrete struc- 
tures subject to high temperatures. However, much less attention has been given to the qualitative 
properties of the model. 

Let us briefly present the single-fluid-phase, purely macroscopic, model for prediction of hygro- 
thermal behavior of heated concrete. The one-dimensional heat and mass transport in concrete 
wall is governed by the following parabolic system [H [51 [7] : 
energy conservation equation for concrete: 

{a,we + psCse)t = (A(0, p)e,), + {c,,K{e, p)ep,)^ ~ h^du (i) 

water content conservation equation: 

vut = {K{e,p)p,), + dt. (2) 

The primary unknowns in the model are temperature 9, water content w (the mass of all free 
evaporable water per m'^ of concrete) and pore pressure P. In order to complete the system ([T])-®, 
let us note that the water content w is connected with temperature T and pore pressure P via 
sorption isotherms $ = P/Psat = ^{9,w) (<& denotes the relative humidity and Psat is saturation 
vapor pressure), which need to be determined experimentally. The thermal conductivity A and 



permeability k are assumed to be positive smooth functions of their arguments. Further, Cw 
represents the isobaric heat capacity of moisture, ps and Cg, respectively, are the mass density and 
the isobaric heat capacity of solid microstructure (excluding hydrate water), d is the total mass of 
the free water released in the pores by dehydration and hd denotes the enthalpy of dehydration per 
unit mass. 

To complete the model, one should prescribe the appropriate boundary conditions across the 
boundary and initial conditions on primary unknowns. In the case of heat transfer through the 
boundary, the Neumann type or radiative (for heat flux) boundary conditions, respectively, are 
frequently under consideration in practice: 

-Aie, P)e^-n = ac{9-9oo), (3) 
-A(0,F)0, -n = ac{9 - e^) + e<T{e^ - et). (4) 

respectively, in which ac designates the film coefficient for the heat transfer, and doo is temperature 
of the environment. The last expression in Q expresses the radiative contribution to the heat 
flux, quantified by the Stefan-Boltzmann law in terms of the relative surface emissivity e and the 
Stefan-Boltzmann constant cr and the temperature difference {9'^ — 0^). 
The humidity flux across the boundary is quantifled by the Newton law: 

-K{0,P)P,-n = l3,{P-Poo), (5) 

where the right hand side represents the humidity dissipated into the surrounding medium with the 
water vapor pressure in the environment Poo and /3c represents the surface emissivity of water. 

The system accompanied by the appropriate boundary and initial conditions, represents 

a challenging mathematical problem. Models of hygro-thermal transport processes in concrete are 
associated to systems of doubly and strongly nonlinear parabolic equations of type (written in 
operator form) 

dtB{U) - V • A{U, VU) = T{U, VC/). 

Although there is no general existence and regularity theory for such problems, some partial out- 
comes assuming special structure of operators A and B and growth conditions on T can be found 
in the literature. The existence of weak solutions subject to mixed boundary conditions with homo- 
geneous Neumann boundary conditions has been shown by Alt & Luckhaus in [3]. They obtained 
an existence result assuming the operator B in the parabolic part to be only (weak) monotone and 
subgradient. This result has been extended e.g. by Filo & Kacur in [TS], who proved the local 
existence of the weak solution for the system with nonlinear Neumann boundary conditions and 
under more general growth conditions on nonlinearities in U. These results, however, are not appli- 
cable if B does not take the subgradient form, which is typical of coupled heat and mass transport 
models. Thus, the analysis needs to exploit the specific structure of the problem. In this context, 
Vala in [33 proved the existence of solution to the purely diffusive hygro-thermal model with non- 
symmetry in the parabolic part but with unrealistic symmetry in the elliptic term. In 'ST and '25J , 
Li et al. study a coupled model for heat and mass transport arising from textile industry, which is 
described by a degenerate and strongly coupled parabolic system. They prove the global existence 
for one-dimensional problem using the Leray-Schauder fixed point theorem. 

Additional results are available for the time discrete setting. Dalfk et al. [T3] analyzed the 
numerical solution of the Kiessl model for moisture and heat transfer in porous materials. They 
proved some existence and regularity results and suggested an efficient numerical approach to the 
solution of the resulting system of highly non-liner equations. However, the Kiessl model is valid for 
limited temperature range only and as such it is inappropriate for high-temperature applications. 
In 0, Benes et al. extended the work [13^ by proving the existence of an approximate solution for 
the Bazant-Thonguthai model, arising from the semi-implicit discretization in time. 

In the present paper we prove a global-in-time existence of a weak solution to a fully nonlinear 
coupled parabolic system with nonlinear boundary conditions modeling heat and moisture transport 
in concrete walls at high temperatures based on the single-fluid-phase model introduced by Bazant 
& Thonguthai [5]. The main result is proved by means of a fixed point argument based on the 
Leray-Schauder approach and approximation procedure and then carrying out the passage to the 
limit. 

Outline of the paper. The paper is organized in the following manner. In Section[2l we introduce 
basic notation and the appropriate function spaces, present the classical formulation of the problem 
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under consideration and specify our assumptions on data and structure conditions in the model. In 
Section[21 we formulate the problem in the variational sense and state the main result of the paper - 
the global-in-time existence of the weak solution. In SectionHJ we derive the a priori estimates for an 
approximate solution of the auxiliary regularized problem. The solution of the regularized problem 
is obtained by the Leray-Schauder fixed point theorem. Using the limiting procedure we prove the 
existence of the weak solution to the original problem. In Section [SI numerical experiments are 
performed to present the moisture migration, temperature distribution and pore pressure build up 
in the model of concrete wall one-side-exposed to various fire scenarios. 



2 Preliminaries 

2.1 Notation and some function spaces 

Vectors, vector functions and operators acting on vector functions are denoted by boldface letters. 
Throughout the paper, we will always use positive constants c, ci, C2, . . . , which are not specified 
and which may differ from line to line. As usual, for a function (/) = (f>{x,t), (j)x and (j)t indicate 
the partial derivatives with respect to spatial variable x and temporal variable t. Let T > and 
£ > be the fixed values, Vl = (0,£), / = (0,T), Qt = x /, 17 x (0,i] {t > 0). We denote 

by W''P = / > and 1 < p < oo and, especially, = W"'P, where W'''P{n), fc > 0, 

p G [1,-|-(X)], denotes the usual Sobolev space with the norm || • ||vi/''.p(n) and, in addition, LP{fl) 
denotes the usual Lebesgue space equipped with the norm || • ||Lp(n)- 

Let X be an arbitrary Banach space with the norm || • ||a' {X* represents the dual space to X). 
Let r G [l,oo). As usual U{I\X) and L°°{I;X) denote the Banach spaces 

|(/); 4>{t) e X for almost every t e I, J dt < ooj 

\ (j); 4){t) e X for almost every i G /, ess sup ||0(t)||x < oo > 

iI'^iil.(7;x) --^ {f wmwx dt) 



and 



with the norms 



l/r 



and 



II0IIl°°(/;X) := esssup ||0(t)||x- 



Further we introduce the space C(/; X) - space of functions cj) '■ [0, T] — > X, continuous, for which 

\\<P\\c{i;X) ess sup \\(l){t)\\x- 

t6[0,T] 

Define the spaces 

V^^\Qt) = {0 e i'(QT)| e L\Qt)] and V^^I(Qt) = V^'^Qrf ■ 

We will often use the following well-known embeddings that are consequences of Aubin-Lions lemma 
and interpolation inequalities: 

VI\Qt)-^-^L\I-W'-^), (6) 
L\I- W^^\n)) n L°°(/; L\n)) ^ L^{Qt), (7) 
V^^\Qt) C{Qt). (8) 

2.2 Classical formulation of the problem 

Incorporating the relation P = Psat{0)^{d ,w) = V{9,w) (via sorption isotherms $ = P/Psat = 
^{6,w)) into the system ((T|)-([S]) we can eliminate the unknown field P and consider the problem 
with only two unknowns 9 and w. Consequently, the classical formulation of the problem we are 
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going to study reads as follows: 



wt - dt 

{C^w9 + psCs9)t + hddt = {X{0, w 

1 



dt = 



id-de,i9)) 



[6^,(9, w)wj: + Se{9,w)9x] 
[X{9,w)9,] 
-5w{9, w)wx - Sg{9, w)9x 
[-X{9,w)9^[ 



x=0 






/3c {ri9{i,t),w{i,t))-Po,) 
{ac + e(7\9{e,t)\^)9{£,t) - ^{t) 



d{x,0) = 
e{x,0)^9o{x) 
w{x, 0) = wo{x) 



(9) 
(10) 

(11) 
(12) 

(13) 

(14) 

(15) 

nfl, (16) 
nn, (17) 
nfl. (18) 



n Qt, 
n Qt, 

n Qt, 

nl, 

n/, 

n/, 

n/, 



Here we assume that all functions are smooth enough. The unknowns in the proposed model are 
temperature 9, water content w and the function d (the total mass of the free water released in the 
pores by dehydration). A, 6g and are diffusion coefficient functions depending non-linearly on 9 
and w. In (fTTj) r represents the characteristic time of dehydration and the function deq = deq{9) 
defines the water mass created by dehydration at equilibrium at the given temperature 9 |12j . In 
([T5|) the function expresses thermal loading given by (cf. (HI) 

m-^aceoo{t)+ea9t,{t). (19) 

Finally, 9q and wq represent the initial distributions of the primary unknowns 9 and w, respectively. 

2.3 Structural conditions and assumptions on physical parameters 

Here we specify assumptions on material coefficients and data in the model ([^- p^ . 

Ai We assume that the material parameters ps, Cs, C^, hd, r, ac, Pc, Poo, cr and e are real 
positive constants. 

^2 Se, Sw and A are functions, 5^, (Jm, A : ^> M, satisfying for certain positive constants Sq, 
^1, S2, Aq, Ai and A2 and for all z — [zi, 22] £ 



< Sq < 6n,{z) <Si< +00, 

\dz^Syj{z)\, \dzjw{z)\ <S2 < +00, 
< 6q < Se{z) <Si< +00, 

MZ)\<C\Z2\, 

\dM^)\,\dzM^)\ < S2 < 

< Ao < A(2;) < Ai < +00, 

\d,,x{z)\,\d,,\{z)\ < A2 < +00. 



(20) 
(21) 
(22) 
(23) 
(24) 
(25) 
(26) 



^3 deq 



is a function of class and there exist positive constants c?i and d2 such that 



for every z G . 



< deq{z) < dl < +CX), 
Kqiz)\ <d2< +00. 



(27) 
(28) 



Aa We assume V 



for every [^1,^2] G K^- 
A^ is a positive continuous function in [0, T]. 
Aq 9q and wq are positive continuous functions in Q. 



to be a Lipschitz continuous function in such that 7^(1^1, ^2) ■£,2 > 
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3 Main result 

The aim of this paper is to prove the existence of a weak solution to the problem described by the 
system ([ ^ - (jl8p . Now we formulate our problem in the variational sense. 

Definition 1 The triplet [0, tu, d] G L^(/;W^'^) x C (I ; L°° (fl)) is called the weak solution to the 
system ©-(HH]) iff 

{w -d)(j)t + {CwOw + psCsO + hdd) ^ndxAt 

Qt 

+ f {6^{e, w)wx + Se{e, w)ex) </),dxdi + / /3c {v{e{i, t), w{i, t)) - p„o) Hi, t) dt 



+ 1 \{e,w)ex^xdxdt- e{5n,{9,w)wx+5e{9,w)ex)iPxdxdt 

JQt 

i-T 

{{a, + eaW, t)f)e{i, t) - d{t) + c^e{i, t)p, {v{e{i, t), w{i, t)) - p^)] ^(^, t)dt 

= I {wo-do)(t){x,0) + {C^eowo + psCsOo + hddo)ij{x,0)dx (29) 

JO 

holds for all test functions [0, "0] G C°°{Qrp), 4>{x,T) = ip{x,T) = Vx G 51. Here d satisfies the 
equation 

1 j-t 

d{x,t) = - e^'-''>/''deq{9{x,s))ds (30) 

T Jo 

for every t G (0,r) and almost every x G il. 

Remark 2 The coupled problem described by the system (p])- (jl8p . being intensively used in practise, 
does not seem analytically studied so far. To the best of our knowledge, there is no existence result 
for the presented model available. 

The main result of this paper reads: 

Theorem 3 (Main result) Let the assumptions Ai-Aq be satisfied. Then there exists at least one 
weak solution [0,w,d] G L^(/; W^^^) x C(/; L°°(f2)) of the system (in|- (IT51) in the sense of Definition 

m 

4 Proof of the main result 
4.1 The regularized problem 

In the proof of the main result we use the method of mollification in order to get the approximate 
solution of the problem (P))- p^ . Let us introduce the mollified system 

wt-{5^{e',w')wx)x-{Se{e',w)S%)x=dt inQr, (31) 

{c^we + psCse)t + hddt - {\{e' ,w')ex)x 

+ C^{9{6U9',w')wx+Se{d',w)9'x))^ in Qt, (32) 
dt = --{d-deg{e)) in Qt. (33) 

T 

Here 9'^ := * £{9) is a regularization of 6*, £ is a small positive real number, denotes the 
standard moUifier defined for the function 9 (see [1, Paragraph 2.28 and Theorem 2.29]) and £{9) 
is the extension operator extending 9 to be zero on \ Qt (see 1, Paragraph 5.17]). Further we 
set := Ji^^w and in a similar way we introduce 9q, Wq and t?^, respectively, as regularizations of 
6*0, Wq and i?, respectively. 

Substituting ((3T|) into ((32|) . equations ((3T|) - (l33l) can be rewritten as 

wt-{SUd',wnwcc)x-iSe{9',w)9l)^^dt in Qt, (34) 

[Cn^w + psCs)9t + {C^9 + hd)dt = ,w')9x)x 

+ C^{5n,{9',w')wx+5e{9',w)9%)9x in Qt, (35) 

dt^--{d-d,q{9)) inQT. (36) 

r 
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The governing equations are supplemented by the foUowing regularized boundary and initial con- 
ditions 



x=0 



x=0 







in /, 
in /, 
in /, 



in /, 

in i7, 
in 17, 
in fl. 



(37) 

(38) 

(39) 

(40) 
(41) 



[-X{9',w')9x] = (ae + ea\9'{e, t)\'')9{£, t) - r (t) 

X—l 

d(x,0) = 
9{x,Q)^9l{x) 
w{x, 0) = Wq{x) 

In this subsection we prove the existence of a strong solution to the regularized problem 



(42) 
(43) 



Theorem 4 The system p4 p - (|43p has a strong solution [9^,w^] £ Vj' (Qt), such that 

\\[0e,We]\\L2{I;Wi.^) < C, 
\\[0e,Ws]\\L--{I;-L^) < C, 

where the constant c does not depend on e. 



(44) 
(45) 



Proof 4.1.0.1 Proof of Theorem^ follows from the a priori estimates in the "weak" classes of 
functions by the usual Leray-Schauder fixed point arguments. 

Solution to an auxiliary problem. For any given couple [9,w] e L'^{I;W-^''^) and < C ^ 1 
consider the initial-boundary value problem 



Wt - {dU9',w')Wx)^ - {dgi9',wWx^^ = dt 

[C^w + psCs)9t + {C.,,9 + hd)dt - {X{9' ,w')9x)x 



in Qt, (46) 









+ [SU0',w')wx+Sg{9',w)9lj 


9x in Qt, 


(47) 






dt 


^~-{d~Cd,,{9)) 

T 


in Qt, 


(48) 


[5^{9' ,w')wx + 5g{9' ,w] 




x=0 


= 


in /, 


(49) 


[A(r>^) 


9x] 


x=0 


= 


in /, 


(50) 


-5w{9'' ,w'')wx ~ 5e{9^,w 




x=e 


= I3,{r{9%e,t),w{i,t)) ~ CPoo) 


in /, 


(51) 


[~X{9',w' 


9x] 


x=e 


= ia, + ea\9'{£,t)\'')9{e,t)-C^^t) 


in /, 


(52) 




d{x,0) 
9{x,0) 
w{x, 0) 


= 

= Cwoix) 


in i7, 
in i7, 
in D,. 


(53) 
(54) 
(55) 



The proof of the existence of the solution to the system (|46 p - (|55l) is split into three steps: 

Step 1 First we treat the initial problem (|48p and (|53p. Let d be a solution of the ordinary 
differential equation 

dt + -d= -Cdeq{9) 

T T 

(which holds for almost every t (0,T) and x eQ) with the initial condition 



dix,0) = mQ. 
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Then 

d{x,t) = - I e^'-*'>/^Cdeq{d{x,s))ds (56) 



for every t £ (OtT) and almost every x £ fl. Hence, for almost every x £ n we have d{x,-) £ 
Moreover, d and dt are nonnegative almost everywhere in Qt (d{x, •) is nondecreasing 
with respect to t). 

Step 2 The system of equations (|46p . (|49p . (j5ip and (|55p (written in non- divergence form) is a 
special case of the problem (cf. 123. Chapter V., (7.1)-(7.2)J) 

wt - a{x,t)wxx -h{x,t,w,Wx) =0 in Qt, 

a(0,t)w^(0,i) + gi(0,t,u;) = in/, 

-a(£, i) + 52 t , = in /, 

w{x,0) — C,Wq{x) in il. 

/n our case we have 

a{x,t) — 5t^,{9'^ ,w'^), 

h{x,t,w,wx) = (^5^{e',w'))^Wx + (^5e{e'M)jl + 6e{e',w)ei^ + dt, 

<7i(0,t,H = 5e(e'{0,t),w{Q,t))~ei{Q,t), 

52(^, t, u-) = -56{e'{£, t), wi£, t))6l{£, t) - PcV{e%£, t), w{£, t)) + p.^P^. 

Under the assumptions A1-A4 and Aq, following the classical parabolic-equation theory for quasi- 
linear equations 1231, for any given [9, w] G L^{T, W^^^) and d given by (j56p the problem (I46p . (j49p . 
(jSip and (j55p admits the unique solution w G V2'^{Qt) (see '23, Chapter V., Theorem 7.4])- Note 
that V2'^{Qt) ^ C{Qt). Now we prove that w is non-negative in Qt- Let = maxjO, w} and 
w~ = max{0, —w}. Test by w~ to get 



- / |w"(t)pd.-E + / 5^{0'',w^)\w~fdxds+ / dtw~dxds 

+ I PcCPooW-{£,s)ds- f l3cP{0'{£,s),w{e,s))w-{£,s)ds 
Jo 

de{0',w)eiw-dxds. (57) 



Since w = — w we can write 

f r{d'{£,s),w{£,s))w~{£,s)ds = / r{e' {£, s), -w^ {£, s))w- {£, s)ds 
Jo Jo 

= - I V{e'{l,s),-w-{£,s)){-w-{l,s))ds 
Jo 

and due to A4 the last integral on the left-hand side of the equation (j57p is always non-positive. 
Hence, applying (I23p and the Young's inequality to the term on the right hand side and neglecting 
the non-negative terms on the left-hand side, we arrive at the inequality 

1 / \w^{t)fdx+ 1 SniiO" ,w'')\w;^fdxds < [ c\w^ \\ei\\w^ \dxds 

2 Jn JQt JQt 

< c(C)||el||L-(Q,) / \w-\Mxds 
JQt 

+^ll^xlU~(g*) / \w^\^dxds. 



Hence, choosing ^ sufficiently small and taking into account (|20p . we can use the Gronwall lemma 
to conclude w~ = and thus w >0 in Qt- 
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step 3 Now having w G V2'^{Qt), w > 0, we get 9 as the solution of the system (|47p . (j50p . I|52p 
and (|54p . We can write this linear parabolic problem in the form 



9t - ai{x,t)9xx - a2{x,t)9x + az{x,t)9 


= f{x,t) 


in Qt, 


(58) 


9Ao,t) 


= 


in /, 


(59) 


X{9'i£,t),w%£,t))9x{e,t) - {a, + ea\9'{i,t)\^)9{e,t) 




in I, 


(60) 


9{x,0) 




in ri, 


(61) 



where 



X{9'',w^ 



(Cn,w + psCs) 

a2{x,t) = — ^ h{e\w') + G^(5^{9',w')wx + 5e{9',w)~9%) 



a^ix^t) 



{Cn,W + psCs) 

Cwdt 

(CniW + psCs) ' 

-hddt 

{CniW + PsCs) 



and ai £ C{Qrp), 02 G i^+'^(/; L°°(f2)) (e is a small positive number), 03 and f G L°°(il). By the 
linear theory for parabolic problems (see 115, Theorem 2.1]) there exists the uniquely determined 
solution 9 G V^'^{Qt) of the problem ((58l) - (|6T|) . 

Finally, let us conclude that for any given couple [9,'w\ G L^(/;W^'^) and Q < C, < \ we have 
[9,w\ G V2'^((5t) as the solution of the problem (|46p - (j55p . 

Basic a priori estimates. In this paragraph we prove some a priori estimates for 9 and w. 
Test (j34p by Cyj9^ and (j35l) by 29 to obtain (adding these both resulting equations) 

[ {9^{C^w + psCs)) ^dxds + [ 2A(r, u;^)|6'^pda;ds + [ C^9^dsdxds 
■iQt " JQt JQt 

+ / 2hd9dsdxds+ f C^9{£,sfl3,iP{9%£,s),w{£,s))~CPoo)ds 

JQt Jo 

+ f 29{e,s){{ac + ea\9'{e,s)\^)9{e,s)-C'9''{s))ds^0. (62) 
Hence, simple modifications yield 

{x,t)^{C^w{x,t) + psCs)dx + I 2A(6'^^^;^)|6'^pda;ds+ / a,9^dsdxds 

+ [ C^9{£,sfl3cV{9^i,s),w{i,s))ds+ [ 29^{e, s){ac + ea\9%l,s)\^)ds 

JQ Jo 

(x, Of{C^uw{x, 0) + psCs)dx - / 2hd9dsdxds 



+ C C^9{£,sfp,Po^ds + C 29{i,sW{s)ds. (63) 
Jo Jo 

The last three integrals can be estimated using Young's inequality and 118, Remark 4], respectively, 
as follows: 

/ 2hd9dsdxds < hd 9'^dxds + hd dl{x,s)dxds, (64) 
JQt JQt JQt 



and 



C^9{£,sYPcPo.ds < C^^cPoo / e\9,\' + Cie)\9\'dxds (65) 

JQt 

29{£,s)'0'{s)ds<e f 9^ {£ , s)ds + C {e) ( d'{sfds, (66) 
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where e represents sufficiently small positive real number. Note that since w > 0, yields 

c^9{e, sf(3,r{d'{£, s), w{e, s))ds > o. 



Taking into account and neglecting the non-negative terms on the left-hand side in 

(1631) we arrive at the estimate 



ci / 0{x,tfdx + C2 / \e^fdxds + C3 / e^{e,s)ds< Ci + C5 9 dxds. 
Jn JQt Jo JQt 



(67) 



Hence, by means of GronwalVs lemma one checks that 

ll^'llL2(/;wi,2(f^)) < c, (68) 

\\0\\L^(I:mn)) < C (69) 

for the positive constant c being independent of C, and e. 

Now let us derive some uniform estimates for water content w. Multiplying the equation 
by w and integrating over Qt we get 



w{x,OYdx-\- / wdsdxds (70) 

O jQt 



- I w{x,tydx+ / Sn,iO'',w'')\wx\ dxds + / Se^O^ ,w)eiwxdxds 
+ f w{£,s)P,{V{0'{e,s),wi£,s)) - CPoo) ds 



and consequently 



1 / w{x,tfdx+ f 6^(9", w'')\w^\'^ dxds + f w{e,s)pcViO''{e,s),w{£,s))ds 

2 Jn JQt JQ 



< - / w{x,0) dx + — w dxds + — / dg{x,s)dxds 
2 2 Jq^ 2 Jq^ 

6g{e^,w)eiwa,dxds + C f w{e, s) (3cPocds . (71) 
Jo 

By A4 we deduce 

w{e, s)l3cP{9%e, s), w{e, s))ds > 0. (72) 
Further, by (|22p and the Young's inequality we have 

[ Sgie^,w)9lw^dxds < e [ w^dxds + C(e) / \5e{9'' ,w)\^\9l\'^dxds 
JQt JQt JQt 



< e / wldxds + ciC{e)\\9\\l2^j.^i.2^^)), (73) 
J Qt 

while using the Cauchy's inequality and U8\. Remark 4-] one obtains the estimate 

/■* 

w{i,s)l3cP^ds < C1+C2/ w^{l,s)ds 

Jo 

< C1+C2 ewl + C(e)w^dxds. (74) 
jQt 

Now taking (|7ip - (l74p together and using (1201) and (j68p . we arrive at the estimate of the form 

ci / ^(x, t)^dx + C2 / ui^da^ds + C3 / w^dxds < C4 + C5 / w^dxds. (75) 
Jn JQt JQt JQt 

By the Gronwall 's inequality we get the uniform estimates for w 

\\w\\L^{i;W^-Hn)) < c, (76) 

lkl|L~(/;L2(f2)) < C (77) 

for the positive constant c being independent of and e. 

Having established the a priori estimates for the solution of the regularized problem we are ready 
to complete the proof of Theorem using the Leray-Schauder approach. 
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Leray-Schauder fixed point arguments. Denote by X = L^(/; W^'^), take arbitrary [9,w] £ 
X, ^ S [0,1] and define the couple [0,w] G as the solution of the problem (1461) - (|55p . 

Define the nonlinear mapping A : X x [0,1] X, given by the equation [0,w] = A{[0,w,Q). It 
is a technical procedure to check that the mapping A is continuous and compact (cf. (OJ. For 
C = we have A{[9,w,0]) = [0,0] for all [9,w] G X. The estimates (j68p and (j76l) imply that [0,w], 
the solution of the problem [0,w] — A{[9,w,(^]) for some C G [0,1], is uniformly bounded in X. 
Now the existence of at least one fixed point [6^,We] G X, A[9^,w^,l] = [^cWs], follows from the 
Leray-Schauder theorem JWi. Consequently, [0^,We\ = A{[6eTWeA]) G '^'^^iQr) and [O^^w^] is the 
strong solution of the system (l34l) - (j43l) . The estimates (|68)) and ((76)) imply (j44l) . Clearly, (|69)) and 
(1771) yield (j45p . The proof of Theorem is complete. 

4.2 Passage to the limit for e: — i- 

To complete the proof of the main resuh stated in Theorem [3] we pass to the hmit for e — > 
and study the convergence of the solution [6ir,w^] of the system (IM|) - (P51) . We present various 
convergence results based on the uniform estimates (|68| ~ ((69|) and (|76l) -(l77 | and, in addition, prove 
some a priori estimates for the strong solution [0^, w^] of the problem ([M)) - (|151) . which is equivalent 
to the system pip - ([55|) with the boundary and initial conditions ([57 1) - (|^5| . 
First, recall that by Theorem S] we have 

ll^e||L2(/;VFi.2(a)), ||^e||L=°(/:L2(iJ)) < C, (78) 

ll'«£||L2(/;VKi.2(n)), \\we\\L--(i-L^{n)) <c (79) 
and from the equation (j46p we derive that 

\\{We)t\\L'^(I-W^.'^(Q.Y) < C. (80) 

As a consequence of the preceding a priori estimates we see that there exist functions w G L^{I; W^''^{il)) 
and 9 G L^{I]W^''^{fl)), Wt G W^^'^(r2)*), such that, along a selected subsequence, we have 

{sj — s> 0+ as j oo) 

w weakly in £2(7; W^i'2(^^^^ ^g^^ 
{w,^)t Wt weakly in L2(/;H^i,2(^)*)^ (82) 

w almost everywhere in Qt (83) 

and 

61., weakly in L'^{I; W^-'^{n)). (84) 

By the embedding L^{I;W'^''^{n))r\L°°{I; L'^{n)) =-> L^{Qt) we get the uniform bound (using ^ 
and ([751)) 



\\e,{5^{eiwl){w,),+5e{eiw,We),)\\L.n(Q^) 

< Cl||0e||L6(Q,) {\\iw,U\mQr) + II (^s)^ || L^(Q,)) < C2 (85) 

and similarly 

\\CniWe6e + PsCsOe\\L3/2(o^T;W^.^/Hn)) < C- (86) 

From the equation p2p we obtain the uniform estimate 

WiC^w^e^ + psCsOe)t\\Ly^i;W^-^n)') < c. (87) 

Since 

where /3 is a small positive real number, the Aubin-Lions lemma yields the existence of x G 
i^3/2(/; p^i-/3,3/2(^^^ such that (modulo a subsequence) 

C^w.^e,^ + psCs0e, ^ X strongly in 1^/^(1; W'~^^^/^{n)). (88) 

Since ((55)) yields the almost everywhere convergence and w^^ converges almost everywhere to w, we 
conclude 

6^. 9 almost everywhere in Qt- (89) 
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Hence, Cyjiv^jO^. +PsCsdej converges almost everywhere to C^wO + psCsd and x = CwwO+psCsO- 
Finally, dST]) yields 

{C^w,£,^ + psCse,^)t ^ {C^,we + psCse)t weakly in L^'^{I- W^'^'i^)*). (90) 
By dUIl-dHll) and ([89]) we conclude 

\{ei^,wl.){e,.)^^\{e,w)e^ weakly in ^^(Qy)^ (91) 

5^{ei^,wl^){we,)^ + 5g{ei^,w,^){ei^)^ ^ +,5e(^?,?«)0. weakly in L^iQr) (92) 

and finally, 

e,^ (5^{ei^,wl,){we,). + &e{ei^,w,.){ei^)^^ ~^e{5^{e,w)w^+5B{e,w)e,) weakly in L3/2(q^). 

(93) 

Now let us present the convergence of the boundary conditions. By ([75)1 and we deduce 
that there exists a subsequence (not relabeled), such that 

We (0, •) — )■ 111(0, •) and ■) — > w{(.^ ■) strongly in L^{I) and almost everywhere in /. 

(94) 

Taking fixed x = or a; = ^ and using ([55)) with x — C^wO + psCsO we conclude 

C„We^.(0,-)6'e, (0,-) + PsCs6'e, (0,-) ^ C^w(0,-)6'(0,-) + PsCs6'(0,-) almost everywhere in /, 

(95) 

C^w,^{l,-)OeA^r)+PsCsee^{e,-) ^ C^w(A-)0(A-) +psCs^(^,-) almost everywhere in /. 

(96) 

Now imply 

e.^ (0, •) ^ 6'(0, •) and 6,^ {£, •) e{£, •) almost everywhere in /. (97) 

In order to get convergence results applicable to the radiative boundary conditions, we need 
"better" uniform estimates than (I75|) . 

Lemma 5 Let [6g,We] G be the solution of the system p4D - (j43p . Then is uniformly 

bounded in L°°{I;L'^{^)), i.e. 

II^s||l~(/;L4(o)) <c, (98) 
where the constant c does not depend on e. 

Proof 4.2.0.2 Test (|34l) hy C^jO'^ and (|35|) by A9^ to obtain (adding the both resulting equations) 

e{x,tf{C^w{x,t)+ psCs)dx+ ( l2e'^\[e\w^)eldxAs+ I iCy^e^d^dxAs 

jQt JQt 

4M'4dxds+ / C^0(£,s)%(P(r(^,s),u;(£,s)) -CPoo)ds 
t "'0 

Ae{l, sf {{ac + ea\9'^{l, s)f')9{l, s) - Ci9"(s)) ds 
e{x,Q)\C^w{x,Q)+ psCs)dx. (99) 
Simple calculation directly leads to the inequality 

ci / e{x,tfdx + c2 I e^eldxds + c^i f e{e,sfds 

Jn JQt Jo 



< I 0{x,O)'^{C^w{x,O) + psCs)dx - I Ahd0^dsdxds 
n 



+ / C^e{e,sri3^Poods+ / 4e{£,sY'd'{s)ds. (100) 



11 



The integrals on the right hand side of (jlOOp can he estimated using Young's inequality and 
Remark 4], respectively, in the following manner: 



further, 



and finally. 



(101) 



(102) 



/ ihdO^dsdxds <3hd O'^dxds + ha dj{x,s)dxds, 
JQt -/Qt JQt 

^ C^6{i, s)^pcPoods < C^l3,Poo (^j eAO^l + C{t)e^dxd^ 

f A9{£,sfi)'is)ds <e f e^{e,s)ds + C{e) f r(s)''ds, (103) 
Jo Jo Jo 

where e represents sufficiently small positive real number. Now, taking (|100p - (jl03D together, we 
obtain 

ci [ 9{x,t)^dx<C2 + C3 [ e^dxds, (104) 
Jn Jqi 

which yields, applying the Gronwall's inequality, the uniform estimate (j98p . The proof is complete. 
By the Sobolev embedding theorem (see Theorem 4.12] or Theorem 8.1.2]) we have 

w^i'2(^) ^ w^/^^^{n). (105) 

Raising and integrating the interpolation inequahty [1] Theorem 5.2] (e means small positive real 
number) 

||6'e||H'i/4+ = .'i(n) < c||6'e||j^t/'4,4(^)ll^e|lL4(n) ^ (106) 

from to T we get 

• J. \ (l+4e)/6 / V (l+4e)/6 



^ lrellH^i/4+«,4(Q)"'-y ^ II "^111^3/4,4(0) ||i/e 11^4(0) ^ 

< M\a \\Q |l(2-4e)/3 

^ '-lrellL2(/;iy3/i.'i(n))ll'^^llL°=(/;L4(n))- K'^'J') 

Now taking into account ([78| , (l98l) , (jl05l) and (jl07p we arrive at the estimate 

ll^e|lL8/(l+4t)(/.^l/4 + ,.4(f2-)) < C. (108) 

Consequently, we get (along a selected subsequence) 

0^^.(0, ■) -> 61(0, •) and e^^{£, •) ^ e{£, •) weakly in LP(/), 1 < p < 6. (109) 

The strong solution [9^,We] € of the problem (ensured by Theorem|3]) is a 

solution of the system pip - ([55|) with the boundary and initial conditions p7p - (|15|) and satisfies 
the variational problem (corresponding to (|3T|) -(l33 l) and (l37l) - (|43l) l 

- / {We - d^) (j)t + {CwOsWe + PsCsOe + hdde)tptdxdt 

Jqt 

+ f iSUei,wl){w,)^ + 6eieiw,)i0l),)cj), + Xi0lwl){0,),^,dxdt 



Qt 

9, {S^{ei,wl){we)^ + 6e{9l,w,M)^)^, dxdt 

'Qt 

+ [ (3ciVmi,t),w,ii,t)) ~ P^) t) dt 







T 

[{a, + ea\9l{e,t)\'')e,{£,t)~^'{t)] i^{£,t)dt 



+ I a^9e{l,t)l3,{V{9l{£,t),We{£,t)) - P^)i,{i,t)dt 



= / wl(j){x,Q) + {C^9lwl+ psCs9l)i:{x,Q)dx (110) 



J2 
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for all test functions [0, V'] e C°°{Qt), (I>{x,T) = ^p{x,T) = e f]. 

The above established convergences (|ST])-(|H3]), ([TOll-dnil) and (|109p are sufficient for taking the 
limit Ej — > as j oo in piOp (along a selected subsequence) to get the weak solution of the 
system ((9))- p^ satisfying ((29|) . This completes the proof of the main result stated by Theorem [3] 

5 Illustration of transport processes in concrete walls at el- 
evated temperatures 

5.1 Illustrative example 

We assume a concrete wall with a width of 120 mm subjected to fire on one side while on the 
other side, the wall is assumed to be thermal and moisture insulated, see Fig. [T] Three different 
fire scenarios, which represent the time dependency of the ambient temperature on the exposed 
side, are employed: (i) standard fire curve (ISO fire), (ii) hydrocarbon fire curve (HC fire), and 
(iii) parametric fire curve (PM fire), see Fig. [2] The last one is more sophisticated than the others 
since for the PM fire, the temperature is dependent not only on the time of fire but also on the 
specific parameters of a given fire compartment (in our case, the following parameters are assumed: 
qtd = 160 MJm-2, O = 0.12 m^/"^, b = 1000 J m'^ s'^/^ j^-i, fire growth rate: medium, see \Tl\). 
As shown in Fig. [21 this curve also includes a decreasing branch that simulates a cooling phase of 
a fire. More information about the fire scenarios mentioned above can be found in 17J. 



On the exposed side of the analyzed wall, the constant parameters of the boundary conditions 
are set to = 0.019 ms'^. Poo = 1-7542 x 10^ Pa, e = 0.7, a = 5.67 x 10"^ Wm'^K'^, and 
ttc — 25, 50, or 35 Wm^^K^^ for ISO fire, HC fire or PM fire, respectively. The temperature 
6*00 (t) is assumed to be equal to the fire temperature according to Fig. 

The uniform initial conditions 9q — 293.15 K and wq = 71.01 kgm^'^ are considered. 

5.1.1 Material data 

The material properties of concrete at high temperatures are assumed as follows. 

The thermal conductivity of concrete A = A{9,P) is adopted from [THl (46)-(47)], where A^o = 
1.3863 Wm-iR-i, = -0.0007272 K"! (see ^ Tab. 2]). The concrete porosity n = n{e) is 
assumed according to |32j (A7)], with no = 0.1, and the saturation S — S{9,P) according to [Til 
(10)]. The permeability of concrete k = k{9,P) = a{9,P)/g is taken from [5, (12a)-(12b)], where 
ao = 10~^^ ms""'^. The mass of dehydrated water is determined by Eq. (fTT|) . which is adopted 
from [H (M-4)]). Here, the term deq{9) is taken from [2, (1.9)]) with d^^^ i^ k ^ 339 kgm'^ and 
T = 10800 s. The sorption isotherm functions are assumed according to with some modification 
based on [TH (73)] and [16, (29)]. Here, the mass of anhydrous cement per unit volume of concrete 
c = 250 kgm"'^ and the saturation water content at the room temperature wqs — 100 kgm^^. The 
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Figure 1: Analyzed wall. 



Figure 2: Fire curves. 
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parameters h^, ps, Cs and Cw are considered to be constant values, namely hd = 2.5 x 10^ Jkg ^ , 
PS = 2400 kgm-3^ Cs = 900 Jkg-^K-i and = 2080 Jkg-^K-i. 

5.1.2 Numerical procedure 

In order to obtain an approximate solution of the nonlinear model, the well-known Galerkin proce- 
dure can be employed. The spacial discretization is performed by the one-dimensional finite element 
method. We consider linear elements with the element size of 0.0005 m (240 elements in total). 
The time discretization is carried out by a semi-implicit difference scheme. In our case, we assume 
the time step At = 0.5 s. The numerical procedure is described in detail in [S]. 

An algorithm arising from the numerical scheme mentioned above has been included in an 
in-house MATLAB code, which is employed to determine the distribution of temperature, pore 
pressure and water content in the analyzed wall. 

5.1.3 Results 

The results obtained for two different times of fire exposure are given in Fig. [3l It is obvious that 
the distribution of the thermo-hygral quantities in concrete at high temperatures is dependent on 
the type of fire scenario used for simulation. For the nominal fires, where the fire temperature 
is a monotonically increasing function of time (in our case the ISO fire and HC fire curve), the 
temperature distribution across the analyzed wall increases with time for the whole fire exposure. 
Contrary to this for the PM fire, the temperature within the structure follows not only the heating 
phase of a fire but also the cooling period, which may lead to the decrease of temperature in some 
parts of a structure, see Fig.[3]and Fig. 2] This holds also for the peak values of other thermo-hygral 
quantities, see Fig. [3] and Fig. S) 

It is obvious, that the usage of the PM fire curve (or other natural fire scenario) provides a 
better prediction of a real behavior of a structure exposed to fire compared to the results obtained 
by assuming some of the nominal fires (e.g. ISO fire or HC fire curve). On the other hand, the 
nominal fire curves are still widely used due to their simplicity and also due to the fact that, in 
most cases, they lead to conservative results. 

As shown in Figs. |3] and HI the rapid heating of a structure and the related moisture migration 
induce the increase of pore pressure near the heated surface. In some cases, this pore pressure 
build up may lead to the spalling of concrete surface layer, and hence, to the eventual collapse of a 
structure, see [9j. 
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